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Abstract 

We demonstrate that configurational electronic entropy, previously neglected, in ab initio ther- 
modynamics of materials can qualitatively modify the finite-temperature phase stability of mixed- 
valence oxides. While transformations from low-T ordered or immiscible states are almost always 
driven by configurational disorder (i.e. random occupation of lattice sites by multiple species), in 
FeP04-LiFeP04 the formation of a solid solution is almost entirely driven by electronic, rather 
than ionic configurational entropy. We argue that such an electronic entropic mechanism may be 
relevant to most other mixed- valence systems. 

PACS numbers: 64.75. +g, 65.40.Gr, 82.47.Aa, 71.15.Mb 
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First-principles prediction of a crystalline material's phase diagram based on the density 
unctional theory (DFT) is a prime example of the achievement of modern solid state physics 
l|. A pure DFT approach is applicable to zero-temperature (zero-T). To study finite-T phase 
stability, one has to identify carefully all the excitations and degrees of freedom involved 
in creating entropy. Typically in alloy theory the focus is on the configurational disorder 
(substitution of different elements or vacancies (V)) while the electronic degrees of freedom 

fin 

are, in the spirit of the adiabatic approximation, integrated out [1|, |2( (phonon contributions 
may give curative corrections gQQ, especially ia systeras with exotic electron-phonon 
coupling [6j, but they are relatively composition insensitive and will not be discussed here). 
For example, many phase diagrams can be satisfactorily reproduced by considering the 
configurational entropy of two elements [l| or element and vacancy Q. Electronic entropy 



is usually t 



lought of as a small quantitative correction and can be calculated from the band 



structure 



[9J: 

'band 



S: 



-k B j n(f\nf+(l-f)\n(l-f))dE, (1) 

where n and / are the density of states and Fermi distribution function, respectively. Only 
electrons within ~ k B T to the Fermi level participate in the excitations, so S^ and is usually 
small. A different type of electronic entropy could arise if electrons/holes (e/h) are localized 
and contribute to the total entropy in the same fashion as the ordering of atoms. One 
would expect such configurational electronic entropy to be particularly important in mixed- 
valence transition metal oxides. Many technologically important materials, such as doped 
manganites, high-T superconductors, Na- and Li-metal oxides, and mixed conductors, fall 
in this category. Little is known about the contribution of localized e/h to finite-T phase 



stability, though previous evidence exist in doped superconductors [10 



111 ] and perovskites 



121 ] that a configurational electronic entropy term (assuming random e/h distribution): 

^loc.rand = [ x]nx+ (l_ x ) ln (! _ x )] j ( 2 ) 

helps explain the entropy of oxidation/reduction. In Eq. [2] a; is the concentration of localized 
electrons or holes. While S 1 ° c,ra,nd can potentially be as significant as the configurational 
entropy of ions, there exists currently no clear demonstration that electronic entropy can 
qualitatively modify finite-T phase diagram. 

In this letter we investigate the effects of configuration-dependent electronic entropy. We 
go beyond a random model such as Eq. [2] and sample electron configurations explicitly. We 
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FIG. 1: The LiFeP04 structure shown with: a) PO4 (purple) and FeOg (brown) polyhedra as 
well as Li atoms (green) b) adjacent layers on Li and Fe sub-lattices, projected along axis a, with 
nearest-neighbor (NN) inter- and intra-lattice pairs highlighted. 

focus on the Li x .FeP04 system. While its high intrinsic Li + mobility makes it of interest 
as the next-generation cathode for rechargeable Li batteries jl^ . it is also crucial to ensure 
good phase equilibration, even at room temperature (RT). So LLRePCU is a good system to 
benchmark theory against. We show that excellent agreement with the experimental phase 
diagram can only be achieved by taking into account configurational electronic entropy, and 
qualitative discrepancies occur if the electron degree of freedom is ignored. 

LiFeP04 has an olivine-type structure with an orthorhombic unit cell. Li removal at RT 
occurs through a miscibility gap between triphylite (T) LiFePC^ and heterosite (H) FeP04 
[l3l | with both phases having a very limited amount of solubility (vacancies + holes (Fe 3+ ) 



|14j. Recent higher-T investigations of the 



in T and Li + ions + electrons (Fe 2+ ) in LP 
Lia-FePC^ phase diagram by Delacourt et al 15| and by Dodd et al [l6j confirm the low-T 
immiscibility, but also find an unusual eutectoid point at 150°C [ijj or 200°C [3] where the 
solid solution (SS) phase emerges around x ~ 0.45 — 0.65. Above 300-400 °C SS dominates 
all compositions (Fig. Gk)- 

This phase diagram is quite unexpected from a theoretical point of view. First, why does 
the system phase separate at all at low-T? In a simplified picture of a generic oxide Li x MO n , 
the Li + ions repel each other due to electrostatics so that ordered intermediate compounds 



3 



are energetically favorable over phase separation, i.e. segregation of Li + (vacancies) into Li- 
rich (deficient) regions. This is indeed the case in many other materials, in which mobile ions 



181 ] . Secondly, what is the 



and vacancies coexist, e.g. L^CoC^, LLjNiC^, Na^OoC^ 
origin of the complex high-T behavior? Transitions from a two-phase coexistence state to a 
solid solution are typically driven by the configurational entropy of the ions in the SS, with 
a maximum transition T near equiatomic A/B composition. The experimentally established 
phase diagram, shown in Fig. [3^, is unlikely to come from such ionic configurational entropy 
unless the effective Li-V interactions are unusually strongly composition dependent. We 
demonstrate that the topology arises from electron degrees of freedom which stabilizes the 
SS near x ~ 0.5. 



We study the LL 



eP04 phase diagram by Monte Carlo simulations based on a coupled 



201 ]. which is a Hamiltonian that explicitly describes the dependence 



cluster expansion 

of the energy on the arrangement of Li + /V and Fe 2+ /Fe 3+ , i.e. both ionic and electronic 
degrees of freedom. In LizFePC^ the Li + ions and vacant sites sit on an orthorhombic lattice, 
of which one layer is shown in Fig. [lb (large green points). On each side of this Li layer is a 
plane of Fe sites (only one plane shown in small brown points). Representing with Aj = ±1 
occupation of site i by a Li + or vacancy and with e a = ±1 the presence of Fe 2+ (electron) or 
Fe 3+ (hole) on site a, the energy can be expanded without loss of generality in polynomials 
of these occupation variables jiol . I20I : 



-^IA A — J% + Ji\ + JijK^j + Jia\ e a + J ab^a e b + • • • 



(3) 



The expansion coefficients J are called effective cluster interactions (ECI), essentially cou- 
pling constants in a generalized Ising model. In its untruncated form, Eq. [3] is exact and 
includes all multi-body terms within one sub- lattice (Li/V or e/h) and between sub-lattices 
though some truncation takes place in practice. To parameterize Eq. [3] we have performed 



GGA+U calculations for 245 Lij-FePO^ (0 < x < 1 ) configurations with su per -ce 



32 formula units using parameter U — J = 4.3 eV 



2l| and other settings in [21 



22 



Is of up to 
. For each 



i/V configuration, usually more than one e/h configuration were considered. The GGA+U 



22 



23 ] app roach is essential to properly localize electronic states (polarons) in this material 



24j. Removal of the self- interaction through proper treatment of the on-site electron 



correlation of localized <i-electrons in GGA+U has previously shown to accurately repro- 



duce the band gap [25(, lithium insertion voltage 



21 



26. 



271 ] and low-T immiscibility 



22|, 
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FIG. 2: Pair ECI vs. site distance (measured from the sites' ideal coordinates in LiFePO^). The 
circled points correspond to NN Li-Li, e-e and Li-e pairs in Fig. [TJa. 

unlike uncorrected LDA or GGA which incorrectly predict stable intermediate LLjFePO^ 
compounds 22]. Ferromagnetic high-spin Fe ions are assumed. At RT (Li)FeP04 is para- 
magnetic 28| and energetic effects of magnetic ordering are small 22j], so the spin entropy 



(~ A;s[xln5 + {l — x) In 6]) is linear in x at RT, therefore negligible in phase diagram calcu- 
lations. 

Our cluster expansion model consists of 29 distinct ECIs: the constant and the point terms 
with no effect on the phase diagram; 7 small triplet terms, which mainly represent slight 
asymmetry between FeP04 and LiFePO^; and most significantly 20 pair interactions shown 
in Fig. [21 Note that these are effective interactions including the effects of many physical 
factors: electrostatics, screening, relaxation, covalency, etc. The Li-Li ECI (diamond) is 
largest for nearest-neighbor (NN) Li + ions, which repel each other strongly for electrostatic 
reasons. As the pairs are separated further, the repulsion is screened considerably. The 
small negative Ju-u at large distance indicates some mediation of the effective interactions 
by lattice distortions. Roughly the same trend is observed for J e _ e . On the contrary, the 
Li-e inter-lattice ECIs are strong short-range attractions that generally become weaker at 
longer distance. The trend in the three curves is not monotonic, since the ECIs contain 
complex lattice factors beyond isotropic electrostatics. The low-T phase separation can be 
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FIG. 3: Li x FePO"4 phase diagram, a) experimental phase boundary data taken from Delacourt et 
al [3] and from Dodd et al 
c) with explicit Li only. 



16l |; b) calculated with both Li and electron degrees of freedom and 



explained by considering the dominating short-range terms. The Li + ions repel each other 
and so do electrons, while Li-e attractions compete to bind them together: if Li + ions stay 
together then the e~ can bind to more of them. The Li-e attractions prevail partly because 
of the host's geometry: the multiplicity of the NN Li-e ECI, the strongest attraction, is two 
per formula unit, while that of NN Li-Li ECI, the strongest repulsion, is one (see Fig.[T]). We 
therefore conclude that phase separation in Li^FePC^ is mainly driven by Li-e attractions 
in competition with Li-Li and e-e repulsions. This is fundamentally different from a system 
where the electronic mixed valence is delocalized, as in metallic Li z Co02 [7i], thereby making 
the Li-e coupling independent of the Li/V distribution. 
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Monte Carlo (MC) simulations combining canonical (interchanging Li/V or e/h pairs) and 
grand canonical (interchanging Li+e together into V+h and vice versa) steps are carried out 
on a 6 x 12 x 12 super-cell, resulting in the phase diagram in Fig. [3b. Phase boundaries were 
obtained with free energy integration. In excellent agreement with the calculated 

phase diagram features a miscibility gap between FeP04 and LiFeP04, and an unusual 
eutectoid transition to the solid solution phase. The eutectoid temperature is only 20-70 K 
off from Q], and the congruent temperatures are about 100-150 K off. We predict the 
enthalpy of mixing at the eutectoid point to be 8.6 meV/formula unit, consistent with the 



measured lower limit 700 J/mol= 7.3 meV/formula unit for an x = 0.47 sample [161 ]. 

To understand better which physics determines the shape of Fig. [3b, we have also per- 
formed calculations in the more "traditional" way, i.e. to consider only the Li/V ordering 
as the entropy generating mechanism, assuming electrons always occupy the lowest energy 
state for each Li/V configuration. The calculated phase diagram (Fig. [3b) shows a simple 
two-phase region, qualitatively different from experiment but similar to typical immiscible 
systems. The striking difference between Fig. [3b and [3b points to the crucial importance of 
explicitly treating the electron degrees of freedom in excitations and finite- T thermodynamics 
of these mixed-valence systems. 

A deeper analysis of the phase diagram in Fig. [3b requires investigation of the entropy 
driving the phase transition. The total (joint) configurational entropy S(Li, e) of the elec- 
tronic+ionic system can be calculated through free energy integration. To partition the 
entropy into ionic and electronic contributions, we note that 

S(Li, e) = S' (Li) + S'(e) + /(Li, e) , (4) 

where I is the mutual information of the two degrees of freedom, and S'(X) = S(X\Y) = 
Y2 y P(y)S(X\y) is the conditional entropy from the X (Li or electron) degree of freedom, i.e. 
the entropy contribution of X with fixed Y, thermal averaged over the marginal distribution 
P(Y). S'(X) measures how random X can be when Y is fixed. If X and Y are independent, 
5" is exactly the entropy contribution from one degree of freedom. We use S' to compare 
different entropy contributions. In Fig. [4] we show the total and separate entropy along 
the solubility limits of the H and T phases (leftmost and rightmost phase boundaries in 
Fig. [3b, respectively), as well as along x = 0.5 in SS. At low-T the total entropy (bold 
lines in Fig. HbO is small, slightly larger in H than in T. The solid solution phase is far from 
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FIG. 4: Configurational entropy per formula unit, a) total entropy and the sum + S' c ; b) 
separate conditional entropy and S' e . 

random: (1) when it first appears at the eutectoid point, its entropy is a mere 0.3 A;#, (2) 
the total entropy of the H phase exceeds that of SS above about 570 K even though its Li 
content is lower, and (3) up to 900 K, well above the congruent points, the total entropy 
1.1 ks of SS(x = 0.5) is still smaller than (complete random) 25 , g OC,rand (0.5) = 1.39 
The difference between S^Li, e) and S' u + S' e (thin dashed curve of Fig. H^) is the mutual 
information 7(Li,e), indicating how correlated the two degrees of freedom are. Fig. Sb 
shows separate S' u and S' e in dashed and dotted curves, respectively. It is noteworthy that 
in all but the T branches S' e is noticeably larger than S^; S' e dominates the SS phase and 
contributes much more than S^ v At the eutectoid point the mixing entropy driving the 
transition into SS is overwhelmingly electronic: 0.19 ks from S' c vs. 0.05 ks from S' u . A 
qualitative explanation for the larger S' e is that the leading J e _ e terms are weaker than the 
leading Ju-u, and the electron excitation spectrum at a fixed Li configuration is lower in 
energy than the opposite. We therefore conclude, to the extent S' represents a separate 
entropy, that the electron degree of freedom contributes substantially more than Li ions to 
disordering of the system, and that the formation of the solid solution state is driven by e/h 



S 



disorder. To our knowledge, no other examples of electronic entropy- driven solid solution 
have been identified, though electronic entropy driven modification of ordering interactions 
through band entropy has been proposed for Ni 3 V 29j. 

Beyond LiFeP04, our approach and results may help our understanding of other mixed- 
valence transition metal oxides with localized electrons. In oxides both electron localization 
and derealization can occur. For example, a system such as Li^CoC^ is metallic for x < 0.9 



30] and explicit e/h entropy is less crucial. LDA and GGA in which mixed valence states are 
delocalized will be an adequate treatment for such system [7]. On the other hand materials 
in which carriers localize require more careful treatment both for their energy calculation 
(e.g. in GGA+U, SIC methods or DMFT [31]), and for their contribution of the electronic 
degree of freedom to the entropy as demonstrated in the present work. An even more com- 
plicated situation arises in materials where electrons can be exchanged between localized 
and delocalized states, as in Ce [3^]. It should be noted that in our MC simulations, e/h are 
treated as classical particles (but not in the DFT energy calculations). If hopping becomes 
so fast that electron wavefunctions overlap, the notion of localized electrons becomes mean- 
ingless, and it becomes difficult to enumerate the eigenstates over which to sum excitations, 
until one reaches the nearly free-electron limit where the band picture is applicable. It is 
up to further investigation to establish quantitative effects of the localized electron degree 
of freedom in thermodynamics of other transition metal oxides. 
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